Three-dimensional disordered conductors in a strong magnetic field: surface states 

and quantum Hall plateaus. 
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We study localization in layered, three-dimensional conductors in strong magnetic fields. We 
demonstrate the existence of three phases - insulator, metal and quantized Hall conductor - in the 
two-dimensional parameter space obtained by varying the Fermi energy and the interlayer coupling 
strength. Transport in the quantized Hall conductor occurs via extended surface states. These 
surface states constitute a subsystem at a novel critical point, which we describe using a new, 
directed network model. 
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The integer quantum Hall effect is one of the most 
striking phenomena observed in two-dimensional electron 
systems [1]. Central to it is the existence of phases in 
which the Hall conductance is constant over a range of 
values for the Fermi energy. In a quantum Hall phase, 
electron states at the Fermi energy are Anderson local- 
ized within the bulk of a sample, but there exist extended 
states at the edge of a sample, which are robust against 
scattering by disorder. It is natural to ask whether 
quantum Hall phases and edge states are unique to two- 
dimensional electron systems, or whether they have ana- 
logues in three-dimensional conductors. 

An obvious way to approach this question is to consider 
a conductor consisting of layers perpendicular to the ap- 
plied magnetic field, each of which, in isolation, would 
exhibit the integer quantum Hall effect. If the inter-layer 
coupling is weak, it is reasonable to anticipate, with in- 
creasing Fermi energy, the sequence of phases sketched 
in Fig.l: insulator, metal and quantized Hall conductor. 
Our aim in the following paper is to investigate theoret- 
ically this phase diagram and, in particular, the nature 
of surface states in a three-dimensional quantized Hall 
conductor. 

Studies of three dimensional conductors in quantizing 
magnetic fields have an extensive history [2]. A variety 
of situations can be engineered in layered semiconduc- 
tors. Multi-quantum- well structures with thick barriers 
between the wells represent the limiting case of uncoupled 
layers, and simply constitute a number of independent 
two-dimensional systems in parallel [3]. By contrast, the 
influence of interlayer coupling is probed in superlattices 
with appreciable dispersion of the electronic spectrum 
in the direction perpendicular to the layers. Accurately 
quantized Hall plateaus are observed [4], as well as an 
oscillatory variation, with inverse magnetic field, of the 
transverse and longitudinal diagonal elements of the con- 
ductivity tensor, fT:,.,:,., and cTzz, suggesting an alternating 
sequence of quantized Hall phases and metallic phases. 



as in Fig.l. It is possible that some of these features 
persist in homogeneous semiconductors. In narrow gap 
semiconductors in the strong magnetic field limit, a tem- 
perature independent a^y in conjuction with transverse 
and longitudinal conductivities, a^x and a^^, that both 
decrease at low temperature [5,6], has been interpreted 
as an incipient Hall plateau. The foregoing examples 
are of particular interest in the present context, since 
it is likely that they show mainly the effect of disorder 
on single-particle motion. In other settings [2], notably 
in the spin-density wave phases of Bechgaard salts [7], 
many-body correlations play an essential role. 




FIG. 1. Schematic phase diagram and density of states, 
p{E). ill energy. E. for the lowest Landan levels of a layered, 
three-dimensional conductor in a strong magnetic field. 



The theoretical treatment of electrons in a disordered 
potential and a magnetic field has appealing simplifica- 
tions in the adiabatic limit, reached if the potential is 
smooth and the field is strong [8]. These have been ex- 
ploited for three-dimensional conductors by Azbel [9]. 
Classically, there are three components to the dynam- 
ics, with widely-separated time scales. In the adiabatic 
limit, the action associated with each component is in- 
dependently conserved. The shortest timescale is the pe- 
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riod of cyclotron motion about an instantaneous guiding 
center: oscillations of tlie guiding center parallel to the 
magnetic field set the intermediate scale: and at long 
times the trajectory followed in these oscillations drifts 
within the plane perpendicular to the field. Quantiza- 
tion of cyclotron and guiding center oscillations leaves a 
reduced problem, involving only guiding center drift. 

Our starting point is a simplified model for the quan- 
tum mechanics of this guiding center drift in layered, 
three-dimensional conductors. It is a natural generaliza- 
tion of the two-dimensional network model for the quan- 
tum Hall effect [10]: each layer of the conductor is rep- 
resented by a copy of the two-dimensional model, and 
adjacent layers are coupled. Two such coupled layers 
have been studied previously [11], as a representation of 
a spin-degenerate Landau level, but the behavior of many 
layers together has not been investigated before. In de- 
tail, the three-dimensional model consists of a network 
of links, each carrying probability flux in the direction 
of guiding center drift, which meet at nodes, where flux 
is scattered between them. Every link is characterized 
by the phase shift that an electron acquires on transvers- 
ing it, and randomness is introduced by choosing these 
phases independently from a uniform distribution. For 
simplicity, the links are arranged on a regular lattice, as 
illustrated in Fig. 2. Scattering at a node can be specified 
by a transfer matrix which relates ingoing and outcoming 
amplitudes (Ain. Aout) on one side of the node, to those 
(Bin-Bout) on the other. With an appropriate choice of 
gauge [10] 



Ain\^( cosh(6») sinh(e) \ ( Bout 
Aout j \ sinh(6») cosh(6») ) \ Bi^ 



(1) 



We choose the scattering parameter, 6, to be the same 
at all intralayer nodes, with a value. 0i. related to the 
Fermi energy, E, by [12] E = ln(sinh^(^i)). Similarly, we 
take a second common value, 02, at all interlayer nodes, 
the tunneling amplitude being t = tanli(^2)- The model 
therefore has a two-dimensional parameter space: t). 



We have investigated the phase diagram of this model 
numerically, using standard transfer matrix techniques to 
calculate the localization length in quasi-one dimensional 
samples, and a finite-size scaling analysis to extract the 
bulk behavior [13]. We study systems of cross-section 
M X M and length L for M < 12 and L < 8.10^. obtain- 
ing Lyapunov exponents with statistical error < 1.5%. 
To concentrate initially on the properties of bulk states, 
we apply periodic boundary conditions in the directions 
transverse to L. As a simple test of our approach, we 
have checked that we obtain the same phase diagram 
from calculations with the layers of the model arranged 
either parallel or perpendicular to L. 

The results are displayed in Fig. 3. Without interlayer 
coupling (t = 0) we reproduce properties of the two- 
dimensional model: extended states exist only at the 
Landau level center (E = 0). Non-zero coupling (t > 0) 
gives rise to a band of extended states, having a width 
in energy, W(t), that increases with t. For t <C 1, one 
expects [14] W{t) oc t^/'^^'', where i'2d is the critical ex- 
ponent for the divergence of the localization length. ^2(ii 
in an isolated layer, and this form is consistent with our 
data. Bulk states in both the low- and high-energy tails 
of the Landau level remain localized even for the largest 
interlayer coupling investigated. We note that, because 
the network model omits inter-Landau level scattering, it 
cannot capture behaviour in the strong disorder or weak 
magnetic field limits, when extended states presumably 
levitate in energy [15]. 
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FIG. 2. Tlio 3(1 network model. Full linos represent links, 
wliirli carry prohahility flux in tlio direction indicated by the 
arrows. Daslied lines represent nodes. Ingoing (outcoming) 
links on the x-z face are denoted by o (•). 



FIG. 3. Phase diagram of the 3d network model. Data ob- 
tained at the points marked /, M and Q are shown in Fig. 4. 



In order to investigate surface states in the model, 
we compare the spectra of Lyapunov exponents in sam- 
ples with periodic and hard-wall boundary conditions. 
One expects surface states, present only in the cjuan- 
tized Hall conductor with hard- wall boundary conditions, 
to carry current without backscattering in samples with 
large cross-section (M 1). and hence to be asso- 
ciated with vanishing Lyapunov exponents. We study 
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qiiasi-ono-dimonsional samples with (using axes defined 
in Fig. 2) tlieir long side parallel to y. periodic boundary 
conditions in the z-direction, and either periodic or hard- 
wall boundary conditions in the x-direction. With peri- 
odic boundary conditions, the Lyapunov exponent spec- 
trum is gapless in the metal, and has a gap (of size 
for large M , where £b is the bulk localization length) 
in both the insulator and the quantized Hall conductor. 
There are only small changes in the spectra for the in- 
sulator and metal on switching to hard-wall boundary 
conditions, and we attribute these changes to finite-size 
effects. By contrast, in the quantized Hall conductor 
this switch has a dramatic influence on the distribution 
of Lyapunov exponents, shown in Fig. 4. With hard- wall 
boundary conditions, M of the (positive) exponents are 
small, and decrease with increasing M . while the values 
of the others are little altered, indicating that there are 
M surface states. This interpretation is reinforced by 
examining the eigenvectors of the transfer matrix corre- 
sponding to the M smallest Lyapunov exponents (Fig. 4. 
inset): their amplitude is concentrated overwhelmingly 
near the sample surface. 
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FIG. 4. Spectra of positive Lyapunov oxpononts. 
Vn-n = \...M'^/2. at tliroo points in tlio phase diagram, 
marked I (0),M (o) and Q (A) in Fig. 3, for sample size 
M = 10. Inset: Amplitude distribution for the transfer 
matrix eigenvector with the smallest positive Lyapunov ex- 
ponent, in the quantized Hall conductor. Squares represent 
links of the model and are shaded according to the mean 
probability flux, p carried by that link: black, p > 0.04; grey, 
0.04 > p > 2.10^^ white, 2.10^* > p. 

Next, we focus on the subsystem of surface states by 
considering the high energy tail of the Landau level, 
where the bulk localization length within our model is 
very short. Then each layer supports an edge state, which 
is decoupled from localized bulk states but coupled to 
edge states in adjacent layers. Edge states in different 
layers carry probability flux in the same sense, and we can 
represent the surface using the two-dimensional directed 
network model illustrated in Fig. 5, which incorporates 
randomness via link phases, as in the three-dimensional 



model, and has a single parameter, the tunneling am- 
plitude between layers, t = tanh(02)- This model is 
clearly highly anisotropic, and represents an example of 
directed scattering, studied previously in other contexts 
[16]. Consider a flnite sample in the form of a cylinder, 
of circumference C and height L , with its axis parallel 
to the magnetic held. Charge transport around the cir- 
cumference is rather simple. It is characterized by the 
Hall conductance, and the restriction of scattering to the 
forward direction ensures that this is quantized. Trans- 
port in the direction parallel to the magnetic field is more 
subtle. We find that the system is at a critical point, and 
the average conductance, gzz- depends on the aspect ra- 
tio, L/C. For L/C <^ 1, the surface states have a finite 
conductance per square, an, which is a function of the 
tunneling amplitude, t, and g^z = {C/L)an- Proportion- 
ality of to sample circumference, C, rather than cross- 
sectional area, is, of course, a signature of conduction by 
the surface, rather than the bulk, and is a characteristic 
of the quantized Hall conductor. In the opposite limit, 
L/C ^ 1, the system is quasi-one-dimensional, states 
have a flnite localization length, along the cylinder 
axis, and g^z decreases exponentially with L. The sys- 
tem is revealed to be critical by the fact that ^ (x C, with 
an amplitude ratio A = (,/C. Remarkably, it is possible 
to calculate both dn and A analytically. 



FIG. 5. The directed network model. Full and dashed lines 
represent links and nodes, as in Fig. 2. 

To evaluate do, we express each element, tij, of the 
C X C transmission matrix, t, between ends of the cylin- 
der in terms of a sum over Feynman paths. Because 
of the directed character of the model, self-intersecting 
paths must wind at least once around the cylinder, and 
make no contribution to the sum in the limit L <C C 
Retaining only those paths that do not self- intersect, it 
is straightforward to average over the link phases. 

From this, using the Landauer-Buttiker formula [17], we 
find 

an = {eyh)tyil-t^). (2) 

To obtain the amplitude ratio. A, we first parameter- 
ize the eigenvalues of t^t as cosh~^(//„L/C). with Vi < 
V2 < ... < vc- In the limit L C . the {vn} are 
proportional to Lyapunov exponents, and therefore self- 
averaging; A = . Moreover, one expects [18] for 
n C that f„ = n ■ v\. In the converse limit. C L. 
rigidity in the spectrum of t^t suppresses fluctuations in 
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the and one again expects [18] that Vn = n ■ vi. In 

this case 

(To = lim lim — — >^ cosli^^(Lnz^i/C) = — ■ 

L^oc C-s-oo C h ^ h 

Making tlio conjecture that the value of Ui is the same 
in both limits, we find 

A = ty(i-t^). (3) 

In order checic this result, we have evaluated the am- 
plitude A numerically, using cylinders of size C < 64 and 
L < 10^, and extrapolating ^/C to large C. The coin- 
cidence sliown in Fig. 6 between the numerical data and 
Eq.(3) is striking. 




Iii(t^l-t))' 

FIG. 6. Comparison of tlie results of a numerical calcula- 
tion of the amplitude A (O) with the analytic expression ( — ) 

Finally, we return to the three-dimensional model and 
consider bulk critical phenomena at the transitions from 
the metal to the insulator or quantum Hall conductor. 
Since bulk states are localized in both the latter phases, 
we expect the same critical behavior at each transition; 
within the network model, this is guaranteed by an exact 
symmetry, arising because higher Landau levels are omit- 
ted. We calculate the localization length exponent, v, 
for comparison with earlier results, obtained from other 
models of the metal-insulator transition in three dimen- 
sions and strong magnetic field. We find, using a stan- 
dard analysis [13], v = 1.45 ± 0.25, which is consistent 
with the value ;/ = 1.35±0.15. obtained both for a layered 
system and a tight-binding model. [19,20]. 

In summary, we have introduced a model of a layered 
conductor in a magnetic field and shown by numerical 
simulation that it has the phase diagram of Fig. 1. We 
have identified the surface states of a three-dimensional 
quantized Hall conductor as a critical system, described 
by a second model, the directed network model, for which 
key quantities are calculable analytically. 
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